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We analyze the entanglement spectrum of Laughlin states on the torus and show that it is arranged in towers, 
each of which is generated by modes of two spatially separated chiral edges. This structure is present for all 
torus circumferences, which allows for a microscopic identification of the prominent features of the spectrum 
by perturbing around the thin torus limit. 
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Introduction — The description of condensed matter 
phases using entanglement measures, borrowed from the field 
of quantum information theory, has led to an explosive growth 
of interdisciplinary work [1]. Despite all this interest, there 
are few cases where entanglement concepts provide physi- 
cal information that is not obtainable through more conven- 
tional quantities, such as correlation functions. One example 
involves topologically ordered states, for which bipartite en- 
tanglement measures have been shown to be useful probes [2- 
4]. Fractional quantum Hall (FQH) states of two-dimensional 
electrons in a magnetic field stand out as experimentally real- 
ized topologically ordered phases, and have recently received 
renewed intense attention partly due to quantum computation 
proposals based on their topological properties [5]. An in- 
triguing feature of FQH states is that their edges have gapless 
modes, described by chiral Luttinger liquids [6, 7]. In this 
work we study the interplay of two edges, through the study 
of entanglement spectra. 

We focus on bipartite entanglement between two parts (A 
and B) of the system. The entanglement spectrum (ES), {^}, 
is defined in terms of the Schmidt decomposition 

^> = 2e-€«/ 2 |^>®|V?>. 

Here | i/j ) is the ground state, and the states | ipf ) (| ipf )) 
form an orthonormal basis for the subsystem A (B). 

Very recently, ES studies have been used [4, 8] to probe 
edge modes of FQH states. The entanglement between two 
partitions of an edgeless wavefunction may seem at first sight 
unrelated to edge physics. However, studies of ES in non- 
interacting systems [9, 10] have found that the entanglement 
spectrum is also the spectrum of an effective "entanglement 
Hamiltonian" confined to the A region of space, which is not 
identical but similar to the original physical Hamiltonian. If 
this similarity holds for interacting systems, the low-lying ES 
would then show an edge structure, even though the total sys- 
tem has no edge. Refs. [4, 8] analyzed the ES of FQH states on 
the sphere with hemispheric partitioning. The Virasoro multi- 
plet structure of the conformal field theory (CFT) describing 
the edge appear in the low-lying part of the ES. 

In this Letter we present and analyze the entanglement 
spectrum of v = 1/3 Laughlin states on the torus. This choice 




Figure 1 : (Color online) Torus setup for block entanglement compu- 
tations. The lowest Landau level is spanned by orbitals which in Lan- 
dau gauge are centered along the circles shown. The arrows indicate 
the chiralities of the virtual 'edges' created by the block partitioning. 



of geometry gives us access to new physics and new analysis 
tools, compared to the spherical case. The natural partitions of 
the torus are cylinder-like segments with two disjoint edges. 
The ES thus contains the physics of a combination of two sep- 
arate conformal edges. We show that this leads to 'towers' 
in the ES spectrum, when plotted against appropriate quan- 
tum numbers. Even in cases where the two edges have dif- 
ferent spectra, the two spectra combine to form towers. The 
two-edge picture provides significant predictive power, as the 
assignment of only a few edge mode energies enables us to 
construct the remaining ES. 

The torus geometry also allows us to adiabatically connect 
to the "thin torus" limit, which is exactly solvable [11,12] and 
has as ground states the Tao-Thouless crystalline states [13]. 
Many features of the ES can be understood starting from these 
simple states, such as the positions of towers and relationships 
between their energetics. The CFT tower structure persists 
even very close to the thin-torus limit, which by itself is an 
uncorrected product state. 

Geometry and partitioning — We study an TV-electron 
system on a torus with periods Li,L 2 in the x- and y- 
directions, satisfying LiL 2 = 2ttN s (in units of the mag- 
netic length). Here N s = N/v is the number of magnetic 
flux quanta. In the Landau gauge, A = Byii, a basis of 
single particle states in the lowest Landau level can be taken 

/ -1/4 r - 1 / 2 i(P L j+mL 2 )x -(y+ P- j+mL 2 ) 2 /2 

as Ipj = 7T i/4 L 1 1 l^m e 1 e L i J ' ' 

with j = 0, 1, N s — 1. The states ipj are centered along the 
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Figure 2: (Color online) Entanglement spectrum for N s = 36 and L\ — 10. Left panels show the symmetric cut and right panels show 
one of the asymmetric cuts. The origin of Ka is chosen to match the Tao-Thouless state. The blue squares represent numerically obtained 
data. The assigned edge modes are labeled by black dots while the combinations of those edges are marked by red crosses. The script letters 
are microscopic identifiers for the two edges combining to form each tower (see text). The striking correspondence of the red crosses with 
numerical data shows that our algorithm based on the two-edge picture allows the reconstruction of the entire entanglement spectrum using 
only the positions of the black dots. The inset shows a CFT tower formed by two ideal chiral edges, the states labeled with their degeneracies. 



lines y = — 2irj/Li (Fig. 1). Thus the ^/-position is given by 
the x-momentum j. 

Any translation-invariant two-body interaction Hamilto- 
nian, acting in the lowest Landau level, can be written as 

n k>\m\ 

where cj^ is creates an electron in the state ^ m and V& m is the 
amplitude for two particles to hop symmetrically from sep- 
aration k + m to k — m. Hence, the problem of interact- 
ing electrons on a Landau level maps onto a one-dimensional, 
center-of-mass conserving, lattice model with lattice constant 
2tt/Li. The Laughlin states at v = 1/3 are, for all L\, 
the unique zero energy grounds states of the pseudo poten- 
tial interaction = (k 2 - m 2) e -27r 2 (/c 2 +m 2 )/L? [H]. For 
generic interactions, e.g., the Coulomb interaction, the matrix 
elements have a more complicated L\ dependence. Our ES 
data are extracted from ground states of (1), obtained using 
the Lanczos algorithm for numerical diagonalization. 

We bipartition the system into blocks A and B which con- 
sist of I a consecutive orbitals and the remaining N s — I a or- 
bitals, respectively (Fig. 1). Since the orbitals are localized, 
this is a reasonable approximation to spatial partitioning, as 
on the sphere [4, 8, 14, 15]. In this Letter we focus on half- 
partitioning, I a = N s /2 and organize the ES in sectors la- 
beled by the particle number, Na, and the total x-momentum 
(along the block boundary), Ka mod N s , in the A block. 

The v = 1 /3 Laughlin state is three-fold degenerate on the 
torus. The degenerate states are related by translation, and 



correspond to three different thin torus configurations: 

010|010010|010 
100|1 1 0|100 
00l|0 1 l|001 

These states are, for generic (including Coulomb and pseudo 
potential) interactions, adiabatically connected to the bulk 
ground states without gap closing for any L\ [11, 12]. If the 
A partition is taken to be the six middle orbitals, the block 
boundaries are different for the three cases: 0-0 cuts in the 
first case ('symmetric cut'), a 1-0 and a 0-1 cut in the other two 
('asymmetric cuts'). Thus the degenerate ground state wave- 
functions have different ES's (which can equivalently be ob- 
tained by using a single wavefunction, and placing the bound- 
aries at three inequivalent positions). 

Tower structure and CFT identification — Numerical ES 
are shown in Fig. 2 for a 12-particle Laughlin state (N s = 
36). A prominent feature is that the ES consists of 'towers'. 
Most of the towers are symmetric, while some are skewed. 
We interpret each tower as a combination of chiral modes of 
two edges (two block boundaries). An ad hoc assignment of 
a small number of (Virasoro) energies provides the necessary 
input for constructing each tower. 

The number of independent modes of a chiral U(l) CFT 
at momentum k is given by the partition function p(k) = 
1, 1, 2, 3, 5, 7, 11, ... for k = 0,1,2,.... When two linearly dis- 
persing chiral modes combine, one expects an ideal tower of 
states like the one shown in Fig. 2 (inset). The diagonal se- 
quences along the left and right sides are from the individ- 
ual edges, and so have degeneracies p(\ AKa\) at momentum 
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Figure 3: (Color online) (a) The chiral edge levels (type A) identified 
from ES, as a function of torus thickness L\. The rectangle contains 
the single-edge ES levels in spherical geometry [15], here scaled and 
shifted for best comparison with the torus results around L\ ~ 14. 
(b) 'Aspect ratio' of the diamond formed by the four lowest ES levels 
of the A- A tower. The value 2 means a perfect diamond shape. 



shifted by AKa from the tower center. A right-moving mode 
at energy E\ and momentum k\ and a left-moving mode at 
E2 and momentum — k<i will combine to give a mode at en- 
ergy Ei + E2 and momentum ki — k2. The rest of the tower is 
obtained from the states of the two edges through such com- 
binations. Here energies are measured with respect to the vac- 
uum state at AKa = 0. 

The observed towers in the numerical ES can be explained 
by postulating the individual edge spectra to have split de- 
generacies while preserving the Virasoro counting. Two such 
modified edge spectra can be combined to construct towers 
which are less degenerate than the ideal case of Fig. 2 (inset). 
Following this scheme, for the numerical ES towers (Fig. 2) 
we assigned appropriate levels on the right and left of each 
tower to single-edge spectra. Several edges in the Na^N/2 
sectors are identical to edges in the Na = N/2 sector, and 
thus do not need to be independently assigned. As a result, the 
energies of a remarkably small number of single-edge states 
(black dots in Fig. 2) are sufficient to generate the entire ES. 
This is a key result of the present work. 

The assigned single-edge levels have robust relative posi- 
tions for varying torus thickness (Fig. 3a). The relative posi- 
tions correspond well to the single-edge levels extracted from 
ES on a sphere [15], as shown in Fig. 3a. 

Microscopic 'thin-torus' analysis — The adiabatic con- 
nection to the thin-torus (Li — >> 0) limit enables us to under- 
stand features of the ES by perturbing away from this solvable 
limit. In particular, the location and energetics of towers can 
be understood from such microscopic considerations. 

We first consider the symmetric cut and explain the towers 
in the left panels of Fig. 2. The very lowest ES level is found 
at Ka = in the Na = 6 sector (left bottom panel), and 
corresponds to the parent (thin-torus) configuration 

0100100101 010010010010010010 1010010010 . 

1 v __ y 1 

V 

A 

At Li = this amounts to the only entanglement level, and 



for all finite Li its dressed counterpart remains at the very 
bottom of the main tower. The remaining states of this tower 
are all generated from this by Vkm processes (Eq. 1) which 
conserve Na- In particular, the leading levels are generated 
by processes with small k and m. For example, a V21 process 
at the right edge, 0010|0100 0001 1 1000; gives the root 
configuration for the lowest entanglement level at AKa = 1 
and at the left edge this gives the lowest AKa = — 1 level. 
The energetics of the tower is determined by the microscopic 
structure of its edges. We refer to this particular environment, 
10010 1 01001, and the corresponding edge energetics as A. 

Some processes do not conserve Na- For the symmetric cut 
the leading such process is V42 which can change Na by ±1: 
0010|01001 0000|1110. We call such an edge environment 
S, which combines with A type edges to form the observed 
A-B and B-A towers in the Na = 5 sector. Mirror images of 
these exist in the Na — 7 sector. By creating two B edges one 
finds that there are two ways of obtaining the B-B tower in the 
Na = 6 sector, each with momentum shift ±Ia (= ±iV s /2) 
compared to the main tower. The extra two-fold degeneracy 
is seen in our data. The large momentum transfer is because 
a particle leaves block A at one cut, and one enters A at the 
opposite cut. We also predict and observe a nondegenerate 
B-B tower in the Na = 4 sector (not shown). 

We now turn to the asymmetric cut (right panels of Fig. 2). 
Again, the thin torus ground state corresponds to the lowest 
level in the main tower (at Ka = in the Na — 6 sec- 
tor). The two edges, 0100|1001 and 1001|0010, are both de- 
noted A! as they are each others mirror images and hence 
have equivalent energetics. In this case already the lead- 
ing hopping term, V21, changes Na leading to the B' edge, 
0100|1001 0011|0001. Another edge C is generated by 
the F54 process, 1001001|001001 0001111|000001. A 
feature of the ES for the asymmetric cut is that the skewed 
towers within a given N a sector do not in general have mir- 
ror image in that sector. Instead, the mirror images show up 
in the (TV — Na) sector. Note that the energy of the lowest 
AKa = state for each tower (tower vacuum energy) is also 
fixed by the two edges. 

Pursuing the microscopic analysis, one can find many non- 
trivial relations between the energetics at different towers, cuts 
and sectors; we have only outlined the basics. It is also possi- 
ble to derive more quantitative features of the ES as a function 
of Li, through perturbative calculations starting at the thin- 
torus limit. Details will be explored elsewhere. 

Circumference and size dependence — We illustrate the 
roles of N s and Li by focusing on the four lowest levels of the 
A- A tower, which form a diamond shape (e.g., Fig. 4). The 
ratio of the energy of the second AKa = level to the first 
AKa = 1 level, each measured from the lowest AKa = 
level, is an 'aspect ratio' for the diamond shape, and is plotted 
in Fig. 3b. Both the two-edge CFT picture and perturbation 
from the thin-torus limit predict aspect ratio = 2. This is seen 
to hold from the thin-torus limit up to some threshold value 
of Li, which increases as N s is increased. This situation is 
generically true for all CFT features: although the Laughlin 
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Figure 4: (Color online) Comparison of the entanglement spectrum, 
and overlaps, between the Laughlin wave function and the Coulomb 
ground state for L\ — 6, . . . , 14 at fixed N s = 36. Only the central 
part of the most prominent tower for the symmetric cut is displayed. 
We observe that the entanglement spectra of the two wave functions 
become very similar for sufficiently large L\. The appearance of 
"generic levels" beyond the two-edge CFT picture in the Coulomb 
state is indicated by the shaded regions, leading to a tentative notion 
of "entanglement gap" [4]. 

state converges to CFT behavior at any L\ , at larger circumfer- 
ences more particles are required for finite-size convergence. 

For large L\ and finite N s , the edges are close and therefore 
interact, leading to complicated effects such as the aspect ratio 
deviations seen in Fig. 3b. 

Coulomb ground states — The ES for ground states of the 
Coulomb Hamiltonian have more complicated L\ dependence 
(Fig. 4). At smaller L\ only the few lowest levels resemble the 
Laughlin ES; for L\ < 8 the Coulomb ES cannot be generated 
using the two-edge procedure beyond the diamond structure. 
This does however not contradict the larger overlap between 
Laughlin and Coulomb states known at Li — » [11, 16], be- 
cause the higher ES levels are pushed upwards at small L\ 
(c/also Fig. 4). The ES thus exposes correlations much more 
subtle than is visible in overlap considerations. As Li is in- 
creased, more and more Coulomb ES levels match the Laugh- 
lin ES, and the emergent CFT tower structure can be seen. 

Discussion — This work presents entanglement spectra 
calculated through numerical diagonalization, for the Laugh- 
lin state at v = 1 /3 on a torus. We presented two radically dif- 
ferent physical ways of understanding the ES structure. The 
first interpretation is based on a combination of two chiral 
CFT edges. Each of these are individually similar to the edge 
spectrum previously extracted from ES studies on the sphere 
[15]. This interpretation is powerful as it reproduces the en- 
tire ES through the assignment of a few levels. Our second 
approach uses the adiabatic connection to the thin-torus limit, 
and the remarkable fact that the two-edge CFT structure is pre- 
served even close to the thin-torus limit. Perturbative analysis 
based on the simple thin-torus states yields the locations and 
shapes of the towers, and many other quantitative predictions. 

Accessing edge modes in explicit numerical calculations re- 
mains a highly desired but difficult task, due to edge recon- 



struction and other difficulties [7] . Our study of edge combi- 
nations through entanglement calculations provides an alter- 
native track to gaining insight into this issue. 

Our work opens up several important research directions, of 
which we list a few. We expect our results to have interesting 
generalizations to more intricate FQH states, such as the non- 
abelian states. 

Our data at large L\ deviates from the independent-edges 
picture because the edges are close. The present setup thus 
provides the intriguing possibility of studying the interaction 
and interferences between two spatially separated edges, e.g., 
through exploring large L\ features as in Fig. 3b. 

The CFT edge interpretation relies on the idea that the 'en- 
tanglement Hamiltonian' is similar to the physical Hamilto- 
nian. This notion is plausible but entirely unexplored for FQH 
states. There is thus a clear need for constructing and under- 
standing entanglement Hamiltonians. It is also possible that a 
more detailed study of the CFT towers in the ES could yield 
Luttinger liquid features such as the compactification radius 
and more generally the scaling dimensions. 

We acknowledge ZIH TU Dresden and MPG RZ Garching 
for allocation of computing time. 
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